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Abstract 

The magnetic-dipole transition probabilities between the fine-structure levels (ls 2 2s 2 2p) 2 P\/2 — 2 Ps/2 
for B-like ions and (ls 2 2s2p) 3 Pi — 3 P2 for Be-like ions are calculated. The configuration-interaction 
method in the Dirac-Fock-Sturm basis is employed for the evaluation of the interelectronic-interaction cor- 
rection with negative-continuum spectrum being taken into account. The 1/Z interelectronic-interaction 
contribution is derived within a rigorous QED approach employing the two-time Green function method. 
The one-electron QED correction is evaluated within framework of the anomalous magnetic-moment ap- 
proximation. A comparison with the theoretical results of other authors and with available experimental 
data is presented. 

PACS numbers: 32.70Cs 
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I. INTRODUCTION 



ities between the 
Tow- 



During the last years, the precision of measurements of magnetic-dipole (Ml) transitions be- 
tween the fine-structure levels in highly charged ions has been continuously increased 

Since in some cases the Ml transitions are sufficiently sensitive to relativistic - 
correlation and quantum-electrodynamic (QED) effects, this provides good prospects for probing 
their influences on atomic transition probabilities. 

To date, a vast number of theoretical calculations of Ml -transition probabi 
fine-structure levels in highly charged ions has been performed (see, e.g., Refs. 
ever, none of these works have provided a systematic analysis of various effects on the transition 
probability. Such an analysis for the (ls 2 2s 2 2p) 2 P\/2 — 2 P^/2 transition in B-like ions and for the 
(ls 2 2s2p) 3 Pi — 3 P 2 transition in Be-like ions is given in the present paper. 

To calculate the decay rate one requires knowledge of the transition energy and the matrix ele- 
ment of the transition operator. Within this work we employ experimental values of the transition 
energy, which are measured accurately enough for the ions under consideration. 

To analyse the influence of various effects, we decompose the transition probability W l ~*^ into 
several terms, 

W^f = WIT 1 + + + AW£' + AW^l + AWg . 

Here W£?f represents the nonrelativistic Ml -transition probability derived employing the LS- 
coupling scheme. Within the LS-coupling scheme, the amplitude of the magnetic-dipole transition 
is nonzero only between the fine-structure levels and depends on the quantum numbers L, S, and 
J of the initial and the final state J12I1 . This implies that the contribution of the interelectronic- 
interaction vanishes in the nonrelativistic limit. The explicit expression for W£*f is presented in 
Section HH 

The relativistic correction AW^~* is obtained by employing the one-electron Dirac wave func- 
tions for the initial and the final state. For the relativistic case the interelectronic-interaction con- 
tribution is nonzero, but it is generally suppressed by a factor (aZ) 2 /Z. For instance, in case 
of B-like Ar it amounts to about 0.1%. The interelectronic-interaction correction is, however, 
rather important for the (ls 2 2s2p) 3 Pi — 3 P 2 transition in Be-like ions, where the terms 3 Pi and 
1 Pi are strongly mixed. In this investigation two approaches are employed for evaluating the 
interelectronic-interaction correction. The first one is based on the configuration-interaction (CI) 
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method in the Dirac-Fock-Sturm basis, whereas the second one employs perturbation theory with 
respect to 1/Z. Utilizin g th e CI method the relativistic Hamiltonian is specified within the no-pair 
approximation I T3L [14II15I1 . The corresponding contribution to the Ml -transition probability is 
denoted by AW^ . The evaluation of this term is described in Sectionimi 

The no-pair Hamiltonian does not account for the negative-energy excitations in the many- 
electron wave function. However, t : 
basis, can become significant [ 



. ,_>^ L , this effect, being dependent on the choice of the one-electron 

OH. 



In Section lTVi the contribution due to the negative- spectrum 
AW^J is derived. 

In Section|Vj the interelectronic-interaction correction of first order in 1/Z is evaluated within 
a rigorous QED approach employing the two-time Green function method II 1 8TI . Together with ver- 
ifying the terms AW^*^ and AW^J to first order in 1/Z, this provides the contribution AW{^, 
which incorporates the 1/Z interelectronic-interaction corrections of higher orders in aZ. 

Finally, AW^l, is the QED correction. The evaluation of this correction to the lowest orders 
in a and aZ is described in SectionlVD 



The main goal of the present work is to evaluate the lifetimes of the states (2s 2 2p) 2 P^/2 in 
B-like ions and (2s2p) 3 P 2 in Be-like ions to utmost accuracy and to investigate the influence 
of various effects on the Ml -transition probability. The corresponding analysis is presented in 
Section rVTll 

Atomic units (h — e — m — 1) are used throughout the paper. 

II. MAGNETIC-DIPOLE TRANSITION PROBABILITY 

The spontaneous L-pole transition probability from the initial state i to the final state / reads 

Hi 

2tt 



where the initial state has the angular momentum Jj, its z -projection Mj, and the energy E i: , and 
Jf, Mf, Ef denote the corresponding quantum numbers and the energy of the final state. The 
transition amplitude Alm is defined as 



Alm = i L+1 x [^V2L + l(f I Tm I i) . 



(2) 



Here Tfa denote the components of the multipole transition operator T L , which is a spherical 
tensor of rank L. In case of a magnetic transition, T L is proportional to the tensor product of the 



Dirac-matrix vector a. and the spherical tensor C\ f = \/47r/ (2L + 1) Y LM 

Tm = -ij L (ur/c)(a®C L ) L M , (3) 

where ji is the spherical Bessel function and to = E i: — Ef is the frequency of the emitted photon. 
The magnetic transition probability can be expressed in terms of the reduced matrix element of 



1 M 



L 2J 4 + 1 c |U 11 11/1 



For the magnetic-dipole transition (L — 1), the tensor product can be written in terms of the vector 
product 

T = —7=ji wr/c ^ i = — Ji ^r/c M , (5) 

V2 t r 

where /x = — e [r x ot\/2 is the relativistic magnetic moment operator. Taking into account the 
first term in the expansion of ji(ur/c) only and turning into the nonrelativistic limit, one derives 
the following relation between the Ml -transition operator and the magnetic moment operator 

Mnr 

TL = ^-/W- (6) 

3 c 

The nonrelativistic magnetic moment operator is given by 

fx m = -fjiB (L + 2S) , (7) 

where L and S are the orbital and spin angular momentum operators, respectively, and fi B = 
\e\h/2mc denotes Bohr magneton. 

In the LS-coupling scheme, which is realized in the nonrelativistic case, the magnetic-dipole 
transition probability is nonzero only between fine-structure levels with AJ = ±1 0. The 
reduced matrix element of T^ r within the LS-coupling is given by 

(Jf II || Ji) = ~ ~ c ^(Jf II (J + S) || Ji) = - ~l^{Jf II S || Ji) . (8) 



Utilizing the general formula for the reduced matrix element of the spin operator [|2Qfl yields the 
corresponding expression for the transition probability 

4u; 3 

W ™ f = ^Bh i ,L f S Si ,SfSi(S i + l)(2S i + l)(2J f + l){ 7 7 7 } . (9) 
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In particular, for the 2s 2 2p 3 / 2 — > 2s 2 2p 1 / 2 transition one can easily find 

4a; 3 1 

W nT f = gjvl = ^ 2.6973500 • 10 13 [s" 1 ] , (10) 

where A is the transition wavelength, in A. Thus, in the nonrelativistic limit the magnetic-dipole 
transition probability is completely determined by the quantum numbers of the initial and final 
states. 



III. INTERELECTRONIC INTERACTION IN THE BREIT APPROXIMATION 



To evaluate the interelectronic-interaction contributions, we start with the relativistic Hamilto- 
nian in the no-pair approximation, 



= A + HA + , H = J2 h D (*) +J2 V ( 



h3, 



(ID 



i<j 



where h D (i) is the one-particle Dirac Hamiltonian and the index i = 1, . . . , N enumerates the 
electrons. The Coulomb-Breit interaction operator V(i,j) = V c (i,j) + V B (i,j) is specified in 
coordinate space as 

I . , , . .s OLi ■ Otj 1 



V B {i,j) 



(at ■ Vi)(aj • ^jYij . 



(12) 



The frequency-dependent part of the full QED interaction operator, which is beyond the Breit 
approximation and gives rise to the terms of higher orders in aZ, will be considered in Section IV! 
A + is the projector on the positive-energy states, which can be represented as the product of the 
one-electron projectors as 



A + = A+(l) • • • X+{N) 



(13) 



together with 



A+(i) = X;i«n(*))Wi) I • (14) 

n 

Here u n are the positive-energy eigenstates of an effective one-particle Hamiltonian h u 



h Un &n Un 5 



which can be taken to be the Dirac Hamiltonian h D , the 
the Hartree-Fock-Dirac Hamiltonian in an external field 



(15) 

Dirac Hamiltonian in an external field or 
13, 14, 151. 



In order to determine the space of one-electron functions {<f n }ff =1 , we employed the combined 
Dirac-Fock (DF) and the Dirac-Fock-Sturm (DFS) basis set. Here the index n enumerates different 
occupied and vacant one-electron states. For the occupied atomic shells, the orbitals ip n with 
n = 1, . . . , M were obtained by the standard restricted Dirac-Fock (RDF) method, based on a 



numerical solution of the radial RDF equations [21, 22]. Only the Coulomb part Vc(i,j) of the 
Coulomb-Breit interaction operator dT2l was included in the RDF Hamiltonian h DF . 

The vacant orbitals tp n with n = M + 1, . . . , M were obtained by solving the Dirac-Fock- 
Sturm equations 

[h DF -e no ]<Pn = ZnW(r)<p n , (16) 

which can be considered as a generalization of the method proposed in Ref. W23I1 to the relativistic 
Hamiltonian and to an arbitrary constant-sign weight function W(r). For every relativistic quan- 
tum number k we choose an occupied DF function Lp no , which we call as reference DF orbital 
and e no in (fToT) is the energy of this orbital. The parameter £ n in Eq. (IToT) can be considered as an 
eigenvalue of the Sturmian operator. Obviously, for = the Sturmian function coincides with 
the reference DF orbital cp no . If W(r) — ► at r — ► oo, all Sturmian functions (p n have the same 
exponential asymptotics at r — > oo. Therefore, the all set of eigenfunctions of the Dirac-Fock- 
Sturm operator forms a discrete set in the space of one-electron wave functions. The completeness 
of this basis in the nonrelativistic limit is well-known fact. In the relativistic case this problem is 
more complicated and we examined the completeness of the pure DFS basis, which we used in our 
many-electron atomic calculations, numerically, reproducing exact hydrogenlike wave functions 
for the same nuclear charge number Z. It should be noted that the DFS orbitals are orthogonal 
with respect to the weight function W(r) and, therefore, form a linear independent basis set. The 
completeness and linear independence of the combined DF and DFS basis was also examined 
numerically. 

In the nonrelativistic theory the widely used choice of the weight function is W(r) = l/r, 
which leads to the well-known "charge quantization". In the relativistic case, however, this choice 
is not very suitable, since the behaviour of the Sturmian wave functions at the origin differs from 
that of the Dirac-Fock orbitals. In our calculations we employed the following weight function 

= 1 - expHorg] 
{ar) z 

which, unlike l/r, is regular at the origin. 



6 



To generate the one-electron wave functions u n , we used the unrestricted DF (UDF) method in 
the joined DF and DFS basis, 

Un = ^ C mn tp m ■ (18) 
m 

The coefficients C mn were obtained by solving the HFD matrix equations 

FC n = e n SC n , (19) 

where F is the Dirac-Fock matrix in the joined basis of DF and DFS orbitals of a free ion. If nec- 
essary, an arbitrary external field can be included in the F matrix. The matrix S is nonorthogonal, 
since the DFS orbitals are not orthogonal in the usual sense. The negative-energy DFS functions 
were included in the total basis set as well. Eq. ( fT9l was used to generate the whole set of orthog- 
onal one-electron wave functions {u n }^ =1 . 

It should be noted that if even there is no external field in Eq. ( fT9"l) . the set of one-electron 
functions {u n }^f =1 differs from the set of basis functions {<f n }n=v F° r me occupied states, the 
UDF method accounts for core-polarization effects, in contrast to the RDF method. For the vacant 
states the difference is more significant, since the DF and DFS operators are inherently different. 

The many-electron wave function + ( r yJMj) with quantum numbers 7, J, and Mj is expanded 
in terms of a large set of configuration state functions (CSFs) $ a ( JMj) 

V+^JMj) = A + tf( 7 JM/) = J2 c MJMj). (20) 

a 

The standard configuration-interaction Dirac-Fock (CIDF) method is used to find the coefficients 
c a . The CSFs are constructed from the one-electron wave functions u n (fT8t as a linear combination 
of Slater determinants. The set of the CSFs is generated including all single, double, and triple 
excitations into one-electron states of the positive spectrum. 



IV. NEGATIVE-CONTINUUM CONTRIBUTION 

Due to some freedom in the choice of the wave function set {«„}, the positive-energy subspace 
and the corresponding projector A + Eq. (fl4ll can be determined in different ways. This freedom 
can be used to find the optimum many-electron wave function ^ opt within the variational method. 

The energy determined by Hamiltonian (ITTT) can be written as 

E = (* I H np I *) = I H I * + ) , ^ + = A + ^. (21) 
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The real orthogonal transformation (rotation) of the one-electron function space {u n } modifies the 
wave function Q 

^' = exp(T)^ + , (22) 
where the operator T is antihermitian (T^ = — T), 

^ ^ E nm t nm , E nm o) n (x rn c^ rn Qj n . (23) 

n<m 

Here a) n and a n are the creation and annihilation operators of electron in the u n state. The matrix 
elements t nm can be obtained from the variational principle. Then the wave function ^opt satisfies 



the generalized Brillouin theorem H25I1 



(*opt I [4«m, H] | ^ opt ) = . (24) 

This means that the optimum wave function ^ opt is invariable under the single excitations in- 
cluding negative-energy spectrum excitations. However, this does not hold for the wave function 
^> + . Therefore, one should revise the calculation of the matrix element \ A \ of any 
one-electron operator A by admixing the negative-energy spectrum excitations to This is es- 
pecially important for so-called "odd" operators, which mix the large and small components of the 
Dirac wave functions. The Ml -transition operator T 1 © is just of this kind. For this reason, the 
negative-continuum contribution can be significant and depends on the choice of the one-electron 
basis set {u n } 

We consider two equivalent methods for evaluating the negative-continuum contribution to the 
matrix elements of a hermitian one-electron operator A with the wave functions The first 
one is based on the Hellman-Feynman theorem whereas the second one employs the perturbation 
theory. 

The space of the wave functions used to find \& opt is invariant under the transformation U = 
exy (iA), if A is a one-particle operator. Therefore, one can employ the Hellman-Feynman theorem 



1 26] to obtain the expectation value of A 

#(//) = H + fiA, (25) 



— d 

A = o-<*apt(A*) I H{n) | *ppt(/i)) 



At=0 



where it is implied that fj,A is included into the one-particle Hamiltonian, h v (fi) — h u + \xA. Since 
the wave function correction 

5^ = ^ opt - * + = [1 - exp(-T)] ^ opt ~ - E nm t nm ^ opt (26) 

n<m 



accounts for single excitations only, the generalized Brillouin theorem d24b yields 



| H(ji) | * pt(/x)> + (^o P t(/i) I I <$*(//)> = o 



(27) 



and, therefore, 







Neglecting the second quadratic term in the equation above yields 



(28) 



A -TT 



M=0 



(29) 



Thus, the negative-continuum contribution can be evaluated by means of the formula 



AA 



nog 



(tt+(Ai) I # (//) | *+(Ai)> 



A»=0 



(*+| A|*+>. 



(30) 



Alternative expression for this contribution can be obtained employing the perturbation theory. 
Using the equation for the derivative of u n (fi) 





d/j, 



^ (u m (0) | A | u n (0)) 
= 2^ z u m {0) 

IJ=0 -L ^ n ^ m 



we obtain 



AA 



(pos) „(ncg) (u m | A | U n ) 



nog 



2 E E 



(a+a„^ + \H\V + ). 



(3D 



(32) 



n m 

Here the indices (pos) and (neg) indicate that the summation is carried out over the positive- and 
negative-energy spectrum, respectively. 

For the nondiagonal matrix elements, one can derive 





AA^ f 

neg 



and 



64Z = Yr"E 



(pos) ^-^(neg) (u m | A | U n ) 



| a | 



^=0 



(33) 



(34) 



These formulas were used in our calculations of the negative-continuum contribution to the Ml- 
transition amplitude. It was found that the results obtained by means of Eqs. (l33l) and d34l are in 
a perfect agreement with each other. 
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V. HIGHER-ORDER INTERELECTRONIC-INTERACTION CORRECTIONS 



The rigorous QED treatment of the interelectronic-interaction corrections to the transition prob- 
abilities can be carried out utilizing the two-time Green function method [18]. In Ref. H27I1 it was 
done for the 1/Z interelectronic-interaction corrections in He-like ions. Here we perform the 
corresponding calculations for B-like ions. To simplify the derivation of formal expressions, we 
specify the formalis m reg arding the core electrons as belonging to a redefined vacuum (for details 
we refer to Refs. J3Q])- This leads to merging the interelectronic-interaction corrections of or- 
der \/Z with the one-loop radiative corrections. The formulas for these corrections can easily be 
obtained from the corresponding expressions for the one-loop radiative corrections to the transi- 
tion amplitude in a one-electron atom, derived in ifisll . However, the standard electron propagator 
S(e, x, y), which enters the equations, must be replaced by 

S(e, x, y) = S(e, x, y) + 2m £ ^ c (x)^ c (y)5(e - e c ) , (35) 

c 

where the summation runs over all occupied one-electron states refering to the closed shells. Ac- 
cordingly, the total expression is represented by the sum of the pure QED and interelectronic- 
interaction contributions, which correspond to the first and second terms in the right-hand side of 
Eq. (l35l) . As a result, the \ /Z interelectronic-interaction correction to the Ml -transition amplitude 
in a B-like ion between the initial state a and the final state b is 

AA -mt = _ [y^ (fc|/(0)|nc)(n|r^|q) ^ (6|r^|ra)(cn|/(0)|ca) 

1M Vttc ^ I e b -e n e a -e n 

| ^ (bc\I(e a - e b )\an)(n\Tl I \c) (c\Tl(\n)(nb\I(e a - e b )\ca) 

^-^ -I- P _ — P . — P ^— S 





e b + e c - 






(be 


I(£b ~ £ c ) 


cn)(n\Tl I 


a) 




£b ~ 






(be 


I(e a - e c ) 


na){n\Tli 






e b + e c - 







_ ^ \ bc \ 1 y e b ~ £c)\cn)(n\l m\o) _ ^ 



£<j &c &b £n 

(b\Tlj\n)(nc\I(e a - e c )\ca) 



E (bc\I{E a - e c )\na){n\l M \c) _ 
p,_ -L p _ — p — p _ £-~t 



(c|r] f |n)(fen|/(e fc - e c )|ca) 

£<j &c &b &n 



~ \{b\Tli\a) [(bc\I\e b - e c )\cb) + (ac\I'(e a - e c )\ca))\ , (36) 

where 1(e) = a^oFD^e), I'(e) = dl(e)/de, a M = (1, a), and D^ v (e) is the photon propagator. 
In the Feynman gauge it reads 

D^e, x - y) = -4^ J — £2 _ k2 + tQ , (37) 
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where g^ v is the metric tensor. In the Coulomb gauge we have 

Doo(e,x-y) = -: — - — - D iO = D Oi = 0, (i = 1, 2, 3) , 
x — y 

' (38) 
d k exp (zk ■ (x — y)) ( £ kikj 



D ij( e, x - y)= ^J— — +i0 "" [S,„ - -tf j , (i, j = 1,2,3). 

In contrast to Ref. [18], here atomic units and the Gauss charge unit (a = e 2 /hc) are used. 
Expression d36b incorporates the Coulomb-Breit part, which was taken into account by the CI 
method, together with terms of higher order in aZ, the so-called frequency-dependent correction. 
Specifying the operator 1(e) within the Coulomb gauge and setting e = in Eq. (l36l) yields the 
Coulomb-Breit interaction. In this way we can exclude the part, which has already been taken into 
account by the CI method, and obtain the frequency-dependent correction of order l/Z as 

AA iv eq ___ [u_ /o f (bc\AI c (e b - e c )|cn)(w|r^|a) 
1M V 7TC ^— ^ I £b - £ 

(b\T^\n)(nc\AI c (e a - e c )\ca) 



| ^ (bc\AI c (e a - e c )\na){n\T^\c) | ^ (c| | n) (6n | A/ c (g b - e c )\ca) 

n n 

_ ^ (^ c l A/ c(gq - £b)\an) (nlTl^c) _ ( c \ T m \n){nb\ AI c (e a - e h )\ca) 

n n 

+ \{b\T l M \a) [(bc\I' c {e b - e c )\cb) + (ac\I' c (e a - e c )\ca)} J , (39) 

where A/c(ga — £&) = ^c(£a — £&) — ^c(O) and the subscript "C" refers to the Coulomb gauge. 

It should be noted that the total l/Z interelectronic-interaction correction given by equation 
(l36t is gauge independent. This has been confirmed in our calculations to a very high accuracy. 
The calculations were performed employing the B-spline method for the Dirac equation [29]. 



VI. QED CORRECTION 

QED effects modify the transition probability via the matrix element of the transition operator 
and via the transition energy. Since we employ the experimental value for the transition energy, 
we have to consider the QED effect on the transition amplitude only. 

The lowest-order QED correction to the Ml -transition amplitude can be derived by correct- 
ing the operator of the atomic magnetic moment for the anomalous magnetic moment of a free 
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electron. In the nonrelativistic limit it yields 

AV -»• Ma = -A*b [L + 2(1 + « e )S] = /x nr + 5/x a , (40) 

where 

<5/x a = -2/x B K e S , (41) 



0.328 478 965... - + 

27T \7T/ 



(42) 



With the aid of the identity 

(J f || J || Ji) = (J f || (L + S) || J i ) = 5 Jfi j i y/J i (J i + l)(2J i + l), (43) 
one can easily find for the fine-structure level transition (A J = ±1) 

(Jf || Sfl a || Jj) = 2« e (J/ || Mnr || J) • (44) 

Therefore, the QED correction to the Ml -transition probability is given by 

Auj^ 1 

AW ^ = 3^"2jTT (I(J/ 11 (Mnr + W 11 Jl) ' 2 " l(J/ 11 Mnr " Jl)|2) ' (45) 
which yields 

4lj 3 1 

^ 4Ke 3^2jTl l(J/ 11 Mnr 11 Jl) ' 2 " AKeW ™ f ■ (46) 
QED corrections, which are not accounted for by this formula, are suppressed by a small factor 

(aZf. 

VII. RESULTS AND DISCUSSION 

The individual contributions to the Ml -transition probabilities and the corresponding lifetimes 
for B-like and Be-like ions are presented in Tables HI and Hfl respectively. Due to the smallness 
of the E2 transition, which is also allowed, the lifetimes are essentially determined by the Ml 
transition. In case of B-like ions, the experimental values of the transition energy were taken 
from Ref. H for S 11+ , Cl 12+ , K 14 +, Ti 17+ and from Ref. [31] for Ar 13 +. As one can see from 
Table HI the interelectronic-interaction correction AWci turns out to be relatively small due to the 
smallness of the factor (aZ) 2 /Z. The most important contributions are given by the relativistic 
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correction AWn and by the QED correction AWqed- F° r Be-like ions, the transition energies 
were taken from Ref. hi for S 12+ , Cl 13+ , K 15+ , Ti 18+ and from Ref. ||3lt for Ar 14+ . In this case 
the interelectronic-interaction correction A Wei provides an essential contribution to the total value 
of the transition probability. This is due to a strong mixing of the two terms 3 P 1 and 1 P 1 . Except 
for Ar 13+ and Ar 14+ , the uncertainties of the total transition probabilities are mainly determined 
by the experimental uncertainties of the transition energy. For argon ions, the uncertainty comes 
mainly from uncalculated higher-order QED corrections. 

In Table ITITl our results for the lifetime of the (ls 2 2s 2 2p) 2 P^/2 state are compared with other 
calculations and with experiment. It should be noted that the QED correction was taken into 
account in Refs. I10II37I and in the present work only. Besides, different values of the transition 
energy u, indicated in Table ITTTl were used in the different calculations. Since the Ml -transition 
probability W scales as cu 3 , a small deviation in to can change W significantly. For this reason, we 
recalculated the results of Cheng et al. and Froese Fischer lloTl for the (ls 2 2s 2 2p) 2 P 3 / 2 state in 



B-like ions for those transition energies we have employed in our calculations. Table IIVI presents 
these values with (r lloL and without (r° J^]) the anomalous magnetic moment correction and 
the corresponding values (r prcs and Tp res ) obtained in this work. As one can see from the table, 
there is an excellent agreement between our "non-QED" results (Tp res ) and those from Ref. J9] 
(t°). There is also a good agreement between our total results (r pres ) and those from Ref. [10] 
(r). The comparison of our theoretical results with the experimental data shows generally a good 
agreement as well. However, in case of Ar 13+ there is a discrepancy between our 2 P?/2 lifetime 
value 9.538(2) ms and the most accurate experimental value 9.573(4)(5) ms QlSD- 

Table fVl shows a fair agreement of our results for the lifetime of the [ls 2 2s2p) 3 P 2 state in 
Be-like ions with corresponding results obtained by other authors and with experimental data. We 
note that the QED correction has not been considered in the previous calculations cited in the table. 

In conclusion, we have evaluated the magnetic-dipole transition probabilities between the fine- 
structure levels (ls 2 2s 2 2p) 2 Pi/ 2 — 2 Ps/2 for B-like ions and (ls 2 2s2p) 3 Pi — 3 P 2 for Be-like 
ions. The relativistic, interelectronic-interaction, and radiative corrections to the transition prob- 
ability have been considered. Except for a recent high-precision lifetime measurement on Ar 13+ 
iQ,!^] with an accuracy level on the order of 0.1%, most experimental results have large error bars 
greater than 1.5% and, within these error bars, most of them are in a fair agreement with our the- 
oretical predictions. In case of Ar 13+ , the disagreement of our prediction with the high-precision 
experimental value amounts to 0.37% of the total transition probability, less than the value of the 
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corresponding QED correction. At present we have no explanation for this discrepancy. 
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TABLE I: The decay rates W [s _1 ] of the magnetic-dipole transition (ls 2 2s 2 2p) 2 P\/2 — 2 P^/2 an d the 
lifetimes r [ms] of the (ls 2 2s 2 2p) 2 Ps/2 state in B-like ions. Numbers in the parentheses give the estimated 
error. 





oll+ 


/^i 12+ 

CI ^ 


A ..13+ 

Ar JT 


v 14+ 


ll 


Energy [cm x ] 


13135(1) 


17408(20) 


22656.22(1) 


29006(25) 


56243(4) 


w QT 


20.37538 


47.43068 


104.56308 


219.4222 


1599.635 


AW D 


-0.03542 


-0.09302 


-0.23145 


-0.5436 


-5.355 




0.00637 


0.01586 


0.03723 


0.0802 


0.597 


AW ncg 


-0.00159 


-0.00396 


-0.00929 


-0.0206 


-0.176 


AWqed 


0.09451 


0.22001 


0.48502 


1.0178 


7.420 


AW freq 


0.00007 


0.00019 


0.00049 


0.0012 


0.013 


Wtotal 


20.439(5) 


47.57(16) 


104.85(2) 


220.0(6) 


1602.1(5) 


r total 


48.93(1) 


21.02(7) 


9.538(2) 


4.546(12) 


0.6242(2) 
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TABLE II: The decay rates W [s 1 ]of the magnetic-dipole transition (ls 2 2s2p) 3 Pi — 3 P2 and the lifetimes 
r [ms] of the (ls 2 2s2p) 3 P2 state in Be-like ions. Numbers in the parentheses give the estimated error. 





S 12+ 


Q13 + 


Ar 14+ 


K 15+ 


Ti 18+ 


Energy [cm x ] 


9712(14) 


12913(16) 


16819.36(1) 


21571(20) 


42638(4) 


w ni 


12.35488 


29.03947 


64.17056 


135.36899 


1045.4311 


AW D 


-0.02017 


-0.05389 


-0.13242 


-0.31247 


-3.2611 


AW C i 


-0.01302 


-0.04909 


-0.16457 


-0.50484 


-10.0481 


AW nC g 


-0.00053 


-0.00133 


-0.00313 


-0.00704 


-0.0649 


AWqed 


0.05731 


0.13470 


0.29766 


0.62792 


4.8493 


Wtotal 


12.38(5) 


29.07(11) 


64.17(1) 


135.2(4) 


1036.9(4) 


r total 


80.79(33) 


34.40(13) 


15.584(2) 


7.398(22) 


0.9645(4) 
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TABLE III: The lifetimes of the (ls 2 2s 2 2p) 2 P 3/2 level in B-like ions calculated in this work with (r pres ) 
and without (Tp rcs ) the QED correction are compared with previous calculations (rtheor) and experiment 
(r exp ). The lifetime values are given in [ms]. The values of the transition energy [Energy] are presented in 
[cm -1 ]. Numbers in the parentheses give the estimated error. 



Ions 



prcs 



r pros [Energy] 


Tbhcor [Energy] 


Method & Ref. 


T oxp & Ref. 


48.93(1) [13135] 


47.35 [13300] 


MCDF [9] 






49.07 [13115] 


MCBP [101 






49.33 [13144] 


MCDF[33J 






49.07 [13136] 


SS [34] 






49.26 [13122] 


MRCI[35J 






49.60 


RQDO (JJJ 




21.02(7) [17408] 


20.55 [17565] 


MCDF [9] 


21.2(6)[6] 




21.02 [17400] 


MCBP 1101 


21.1(5)[6J 




21.19 [17421] 


MCDF [33J 






21.08 [17410] 


SS [34] 






21.19 [17386] 


MRCI[35J 






21.13 


rqdoujj 




9.538(2) [22656] 


9.407 [22795] 


MCDF [9] 


8.7(511381 




9.515 [22660] 


MCBP [101 


9.12(18)[21 




9.618 [22666] 


MCDF [33] 


9 70(1 5) [41 




9.569 [22653] 


SS [34] 


9.573(4)(5)[81 




9.588 [22657] 


RQDOQJJ 






9.606 [22636] 


MCDF [36] 






9.615 [22619] 


MRCI[35J 






9.534 [22658] 


[37J 




4.546(12) [29006] 


4.509 [29129] 


MCDF [9] 


4.47(10) [51 




4.521 [29044] 


MCBP [Iff 






4.583 [29019] 


MCDF^3J 






4.558 [29004] 


SS [34] 






4.587 [28960] 


MRCIQ5J 






4.577 


RQDOQJJ 




0.6242(2) [56243] 


0.6254 [56275] 


MCDF [9] 


0.627(10)[3J 




0.6150 [56465] 


MCBP [10] 






0.6290 [56258] 


MCDF[33J 






0.6254 [56240] 


SS [34] 






0.6289 [56166] 


MRCIQ5J 






0.6270 


RQDOyjJ 





CI 



12+ 



Ar 



I'M 



K 14+ 



Ti 17+ 



49.16 



21.12 



9.582 



4.567 



0.6271 



MCDF - multiconfiguration Dirac-Fock method 
MCBP - multiconfiguration Breit-Pauli method 
SS - SUPERSTRUCTURE program 

MRCI - multireference relativistic configuration interaction method 
RQDO - relativistic quantum defect orbital method 
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TABLE IV: The lifetimes of the (ls 2 2s 2 2p) 2 P 3/2 level in B-like ions calculated in this work with (r pre s) 
and without (Tp res ) the QED correction are compared with previous theoretical results, recalculated to the 
transition energy (Energy [cm -1 ]) employed in this paper. The lifetime values are given in [ms]. 



Ions 


Energy 


r° 

' pres 


r° (Ref. [9J) 


Tpres 


r (Ref. [40]) 


§11+ 


13135 


49.16 


49.16 


48.93 


48.85 


C1 12+ 


17408 


21.12 


21.11 


21.02 


20.99 


Ar 13 + 


22656 


9.582 


9.581 


9.538 


9.520 


K 14+ 


29006 


4.567 


4.567 


4.546 


4.539 


Ti 17+ 


56243 


0.6271 


0.6265 


0.6242 


0.6223 



TABLE V: The lifetimes of the (ls 2 2s2p) 3 P 2 level in Be4ike ions calculated in this work with (r prcs ) and 
without (Tp res ) the QED correction are compared with previous calculations (r t h CO r) and experiment (r exp ). 
The lifetime values are given in [ms]. The values of the transition energy [Energy] are presented in [cm -1 ]. 
Numbers in the parentheses give the estimated error. 



Ions 



' pres 



rp r es [Energy] ^thcor [Energy] 



Method & Ref. r oxp & Ref. 



?12H 



C1 i3+ 

Ar 14+ 



K 



15+ 



81.16 80.79(33) [9712] 
34.56 34.40(13) [12913] 
15.66 15.584(2) [16819] 



Ti 18 + 0.9689 0.9645(4) [42638] 



83.3 [9743] 


SHF 


80.65 [9720] 


MBPT 


35.7 [12893] 


SHF 


34.60 [12903] 


MBPT 


16.31 [16818] 


MCHF 


16.1 [16824] 


SHF 


15.63 [16834] 


MBPT 


15.76 [16782] 


MCDF 


7.63 [21575] 


SHF 


7.353 [21633] 


MBPT 


0.990 [42653] 


SHF 


0.9615 [42651] 


MBPT 



15.0(7)[1] 
13.4(7)[2J 
15.0(8)[4J 



7.6(5)[5] 



SHF - scaled Hartree-Fock method 
MBPT - many-body perturbation theory 
MCHF - multiconfiguration Hartree-Fock method 
MCDF - multiconfiguration Dirac-Fock method 
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